Ambiguities in the scattering tomography for central potentials 
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Invisibility devices exploit ambiguities in the inverse scattering problem of light in media. Scat- 
tering also serves as an important general tool to infer information about the structure of matter. 
We elucidate the nature of scattering ambiguities that arise in central potentials. We show that 
scattering is a tomographic projection: the integrated scattering angle is a projection of a scat- 
tering function onto the impact parameter. This function depends on the potential, but may be 
multi-valued, allowing for ambiguities where several potentials share the same scattering data. In 
addition, multivalued scattering angles also lead to ambiguities. We apply our theory to show that 
it is in principle possible to construct an invisibility device without infinite phase velocity of light. 

PACS numbers: 42.79.-e, 03.65.Wj, 03.65.Nk, 
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An invisibility device Q, H, IE Q IE IS should guide 
light around an object as if nothing were there. It is 
conceivable that such devices can be made using mod- 
ern metamaterials IE IE ISl- Passive optical devices 
use spatially varying refractive-index profiles for imag- 
ing. Within the validity range of geometrical optics, in- 
dex profiles of isotropic dielectric media are mathemati- 
cally equivalent to potentials for light rays |Ej01- There- 
fore, such an invisibility device corresponds to a poten- 
tial that has the same scattering characteristic as empty 
space. While the inverse scattering problem for waves 
has unique solutions the scattering of rays may be 
ambiguous. Here we show how such ambiguities arise in 
the case of radially symmetric potentials. Our theory 
indicates that it is in principle possible to construct an 
invisibility device where the phase velocity of light does 
not approach infinity, in contrast to all previous propos- 
als for macroscopic cloaking |E IE S - This could inspire 
ideas for developing invisibility devices without anoma- 
lous dispersion p] that could operate in a relatively wide 
frequency window. In addition to applications in a poten- 
tially new area for metamaterials, our theory has wider 
implications for the field of scattering tomography. 

The inversion of the classical scattering in central po- 
tentials is a classic textbook problem that has made it 
into the exercises in Landau's and Lifshitz' Mechanics 
Since Rutherford's experiments, scattering has served as 
an important tool to investigate the structure of mat- 
ter, with modern applications ranging from biomedical 
research to astrophysics. Techniques to infer the struc- 
ture of matter from scattering are often called scattering 
tomography, although, strictly speaking, t hey are not di- 
rectly related to traditional tomography [lOj where the 
shape of a hidden object is reconstructed from projec- 
tions. Here we show that the case of scattering in central 
potentials literally is a tomographic projection in dis- 
guise, but with an interesting twist: the object to be 
reconstructed corresponds to the potential, but may be 
represented by a multi-valued function, allowing for am- 




FIG. 1: Scattering in a ccr 
parameter b is deflected by th( 
U{r) centered at the origin, with r = \J 
and y. The turning point of the trajectc 
Rutherford scattering ir 



al potential. A trajectory i 
.ngle X ill the rotationally 



rith 



+ T/2 in the Cartesian coord 
ory is denoted by rQ . The figu 
1/r potential. 



potential 
dinates x 



biguities. 

Figure n] illustrates the situation typical for scattering 
in central potentials. An incident ray characterized by 
the impact parameter b and the energy E is deflected 
by the angle x- We use polar coordinates with radius 
r and angle (p in the plane orthogonal to the angular- 
momentum vector. The scattering angle is determined 
as [llj 
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Here denotes the turning point of the trajectory given 
by b and E, and p represents the potential as 
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The turning point is given by the largest value of r at 
which the denominator in the integrand of the scat- 
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FIG. 2: RcprcBcntation of the potential by the turning parameter p defined in 
Eq. |2}. The solid line in A shows a potential with a fold in the turning parameter 
where r(p) is multi- valued. Subfigure B shows the eorresponding potential U (r) . 
The dotted lines deseribe p(r) and U{r) for a potential with the same scattering 
eharacteristies. Here r(p) was obtained from Eq. ITol using the definition |^ of 
the seattering funetion W. 

tering angle vanishes, i.e. at which p = b. Reversing 
this relation leads to a physical interpretation for p: p(r) 
describes the impact parameter for which the radius r 
is a turning point. Therefore we may call p turning pa- 
rameter. Figure |21 illustrates the representation of the 
potential using the turning parameter. The largest zero 
of p(r) corresponds to the potential barrier where U = E. 
The potential is repulsive for p < r, zero for p = r and 
attractive for p > r. Note that the inverse function r(p) 
may be multi- valued, as shown in Fig. The addi- 
tional values of r{p) describe the turning points of addi- 
tional bound trajectories for the same energy E and the 
angular momentum that corresponds to the impact pa- 
rameter b. Scattering does not probe such bound states, 
although the trajectories of scattered rays may enter the 
same region for different impact parameters b. As we 
show, the possibility of such elusive bound trajectories 
indicates ambiguities in scattering. 



In the following, we express the description of scatter- 
ing in central potentials as a tomographic projection for 
the integrated scattering angle 
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First, we represent Eq. ^ as 
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in terms of 



W = Hp/r) , a = y/p^ ~ b^ . 
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A prime indicates differentiation with respect to the turn- 
ing parameter. We call W scattering function. When 
r{p) is multi- valued the integration contour is understood 
to follow accordingly. Since 
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we obtain for the integrated scattering angle 
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This result has a simple geometrical meaning illustrated 
in Fig. 13 imagine that a and b constitute a plane of 
impact parameters where one, b, is experimentally acces- 
sible and the other, a, is not. The scattering function 
W depends only on the radius p — y/aP+W, both di- 
rectly by definition |(SJ) and in r{p). Equation Q shows 
that the integrated scattering angle is a projection of the 
rotationally symmetric object W{p) onto the experimen- 
tally accessible impact parameter b in exactly the same 
way as objects are projected in classical tomo graphy |lCl| | 
or Wigner functions in quantum tomography |l2L Il3j | . If 
r{p) is single- valued, one can invert the projection by the 
inverse Abel transformation 
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a special case of the inverse Radon transformation [l^l ■ If 
r{p) is multi- valued one can hide features of the potential 
in the folds of W , as Fig. O illustrates. 

Consider the scattering ambiguities where the scatter- 
ing function W is multi- valued. The simplest case corre- 
sponds to a single fold in W between two turning param- 
eters pi and p2, as shown in Figs. 01 and ^ We use the 
inverse Abel transformation to construct a potential, 
described by W[){p), that exhibits the same scattering 
characteristics as W . Figure |3| indicates that W and Wo 
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FIG. 3: Scattering tomography. The integrated aeattering angle is a pro- 
jection of the scattering function W onto the impact parameter b along the red 
lines of the fictitious parameter a. The radius p in this auxiliary (a, fa) plane is the 
turning parameter The scattering function, defined in Eq. and obtained 

from Fig. |2j, may be multi-valued, as shown here. The folds of W visualize the 
ambiguities of scattering. 
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and p2 where IV_^ denotes the top and W_ the bottom curve of the fold. 

agree for p > pi, because all projections lie under the 
fold. For p < pi the scattering angle x is, according to 
Eq. CEJ, 
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where the integration variable a refers to the turning pa- 
rameter, W follows the solid curve in Fig. 01 with a jump 
at pi, whereas W+ denote the top and W- the bottom 
curve of the fold. Since the inverse Abel transformation 



(|5|) uniquely inverts the first term in x, we obtain for the 
difference between W and Wq 
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by partial integration, utilizing that the boundary term 
vanishes, because W-{p2) = lV+(p2)- Since Wj^ > 
the ambiguous W must exceed Wq in the single- 
valued region inside pi, which implies that the radius 
r = p exp(— W^o) is greater than p exp(— W). The fold of 
multi-valuedness thus magnifies the scattering structure 
of the potential. In particular, for ambiguous scattering 
potentials, the zero of p{r) is closer to the origin than for 
the equivalent non-ambiguous one. Since this zero cor- 
responds to the potential barrier beyond which one can 
hide, nothing is gained, quite the opposite. This feature 
continues in the general case of several folds in W, be- 
cause one could replace W by equivalent single- valued Wq 
with the same scattering characteristics, starting from 
the outmost fold and proceeding to the inside. 

An alternative way of hiding the presence of a poten- 
tial would be to let the trajectories leave at scattering 
angles that are multiples of 27r, i.e to turn them around 
in precisely adjusted loops. Suppose that for impact pa- 
rameters b smaller than a critical bo the trajectories are 
uniformely turned by x = —2711/ and are not affected for 
6 larger than bo. Here v may be a real number, not only 
an integer, for the sake of generality. Assuming that r{p) 
is single- valued, we obtain from the inverse Abel trans- 
formation 

'bo/p+^bl/p^-iy^'^ : p<bo 

■■ p>bo 



(11) 



Figure O illustrates the curves of p{r). Clearly, r{p) is 
single- valued by definition. For < the potential would 
be repulsive, because the trajectories are deflected, but 
in this case the function p{r) itself is multivalued. Con- 
sequently, no central potential exists that uniformly de- 
flects trajectories. For i' > the potential is attractive, 
as one would expect to be necessary for bending trajec- 
tories around the center of force. The case i/ = 1/2 cor- 
responds to a Kepler potential 0| or the Eaton lens 
developed in radar technology. In the limit p we get 
from Eq. (^TJ the asymptotics p/r ~ (2/r)2''/(2''+i) ^ and 
hence, according to Eq. ^ the potential U diverges with 
the power —4:i'/{2v + 1) for small r. One cannot hide 
anything here. In the limit ^ — > oo of infinitely many 
cycles U approaches near the origin the 1 /r'^ potential of 
fatal attraction 

'ill- 

Figure |3|3 illustrates the case where 
the trajectories are turned around by 2tt. 
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FIG. 5: Uniform bending. Subfigurc A shows the turning parameters ob- 
tained from Eq. ITTl for the winding numbers u £ {2.1,0.5,-0.1,-0.2,-0.3} 
and tiQ = 1. For negative u , p{r) is multi-valued and hence unphysical. Subfigure 
B illustrates the uniform loops of trajectories for i> = 1. 

Although one cannot directly apply the ambiguous 
scattering of isotropic and centrally symmetric media to 
construct an invisibility device, one can use their singu- 
larities to improve anisotropic devices. Such a device is 
designed to facilitate a coordinate transformation with 
a hole 0- Anything inside the hole is hidden by con- 
struction 0]. Consider a two-dimensional case in polar 
coordinates. Suppose that the radius r is mapped onto 
r' such that r' reaches the radius of the hole at r = as 

dr' 

ar-" for r - , (12) 

or 

where a and s are non-negative constants. Beyond the 
outer radius 6o of the cloak the coordinates r' shall coin- 
cide with r. Assume in unprimed space the isotropic and 
radially symmetric refractive-index profile n(r) with per- 
fect impedance matching. Reference gives a recipe to 
calculate the dielectric e and magnetic fj, that facilitates 
the coordinate transformation H12|l . We find 

4 = A^r — n-{r) , s' ^ [j! j- n{r) . (13) 

Suppose that we use a profile where corresponds jjj 



to the E—Uoi uniform bending Hll|l with the definition 
(0) and V = 1. If we choose s = 1/3 the singularity of U 
compensates for the zero in the refractive index in real 
space that would otherwise imply that the speed of 
light tends to infinity at the inner surface of the cloak. 
The phase velocity in radial direction is finite. On the 
other hand, the speed of light in angular direction tends 
to zero with the power 4/3. Our simple example indicates 
that invisibility devices with finite phase velocity are pos- 
sible in principle. In our case, wrapping light around the 
invisibility device stratifies the optical wavefronts. How- 
ever, there is a price to pay: light propagation with fi- 
nite phase velocity around an object inevitably causes 
time delays that result in wavefront dislocations at the 
boundary 0. The invisibility is perfect for rays, but not 
for waves. 

Conclusions. — Scattering in central potentials cor- 
responds to a tomographic projection that visualizes 
scattering ambiguities. Such ambiguities are limited, 
though: central potentials are not suitable to achieve the 
same scattering characteristics as empty space. There- 
fore, highly asymmetric refractive-index profiles 0, 0| or 
highly anisotropic media are required to design in- 
visibility devices, which, interestingly, can operate with 
a finite speed of light. Otherwise, trying to hide things 
uniformly from all sides just magnifies them. 
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